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Sacrificial bonds and hidden length in structural molecules account for the greatly increased 
fracture toughness of biological materials compared to synthetic materials without such structural 
features, by providing a molecular-scale mechanism for energy dissipation. One example is in the 
polymeric glue connection between collagen fibrils in animal bone. In this paper, we propose a 
simple kinetic model that describes the breakage of sacrificial bonds and the release of hidden length, 
based on Bell's theory. We postulate a master equation governing the rates of bond breakage and 
formation. This enables us to predict the mechanical behavior of a quasi-one-dimensional ensemble 
of polymers at different stretching rates. We find that both the rupture peak heights and maximum 
stretching distance increase with the stretching rate. In addition, our theory naturally permits the 
possibility of self-healing in such biological structures. 



I. INTRODUCTION 

Many biological, polymeric materials gain their 
strength and toughness through the formation of sacrifi- 
cial bonds and hidden length. Examples include bone 
E], abalone shells [2S| and diatoms [10H13]. Often, 
sacrificial bonds connect two different sites on a molec- 
ular backbone, thereby constraining part of the poly- 
mer from stretching. These bonds are typically weaker 
than the covalent bonds on the molecular backbone; they 
break and release "hidden length" before the molecular 
backbone ruptures. This molecular-scale mechanism has 
been found to greatly increase the total amount of work 
needed to break the material. 

An important example of sacrificial bonds and hidden 
length occurs in the polymeric glue connection between 
collagen fibrils in animal bone [TjHZ] , illustrated schemat- 
ically in Fig. [T] Each intact sacrificial bond shields part 
of the glue strand from contributing to its end-to-end 
distance. Given an end-to-end distance, a glue strand 
of smaller apparent length carries less entropy than one 
with more available length. In other words, the pres- 
ence of sacrificial bonds and hidden length amplifies the 
amount of force that is necessary to stretch the polymers, 
and therefore accounts for the increase in fracture tough- 
ness of the material. Following breakage of one sacrifi- 
cial bond, the corresponding hidden length now unravels, 
causing a force drop as an immediate result of the spike 
in entropy. 

We recently proposed a theoretical model that ac- 
counts for this mechanism and captures the mechanical 
response of the stretched polymer network in the quasi- 
static limit [14]. In that model, we assume that the 
strength of sacrificial bonds is a random variable, pri- 
marily to account for the apparent variability of bond 
strength as observed in several stretching experiments [1]- 
5 . It is well known that the mechanical behavior of 
stretched biological molecules depends on the pulling ve- 



locity. In particular, Bell's theory [TS] implies that the 
maximum force that a molecule can sustain varies as the 
logarithm of the pulling speed, as is observed in experi- 
ments |16| . Our earlier model, however, does not account 
for this rate dependence unless we impose the assump- 
tion that the "random bond strength" distribution itself 
varies as the logarithm of the pulling velocity. Neither 
does our previous model entail a recuperation of strength 
and toughness observed in experiments [3]. 

An understanding of these velocity- and recovery time- 
dependent behaviors is of paramount importance in many 
applications. For example, the propagation of cracks - 
often caused by traumatic injuries in the case of bones - 
is a dynamical process (as opposed to quasi-static). On 
the other hand, self-healing might impede the spread of 
microcracks in bone. To account analytically for these 
behaviors, we borrow the two-state model of protein un- 
folding due to Rief, Fernandez and Gaub |17j . Since 
breakage of sacrificial bonds and protein unfolding both 
involve the forced breakage of noncovalent bonds and the 
unraveling of compact structure, we base our framework 
on the assumption that the kinetics of sacrificial bond 
breakage can be described in a similar manner. The two- 
state model enables us to obtain analytic expressions for 
the transition rate in terms of the force-extension profile. 

The rest of this paper is organized as follows. In 
Sec. II we use the two-state kinetic model to derive condi- 
tions for the breakage and formation of sacrificial bonds 
and show how the pulling velocity relates to these pro- 
cesses. In Sec. Ill we apply the model to a single poly- 
mer chain - a quasi-one-dimensional chain of entangled 
polymer molecules in series. We show that the model 
reproduces the logarithmic dependence of the peak force 
on the pulling velocity and, with a judicious choice of 
the several adjustable parameters, yields force-extension 
profiles that are qualitatively similar to what is observed 
in collagen fibril separation experiments [IHS I16j . In 
Sec. IV we apply the kinetic model to entangled poly- 
mers and examine the effect of the delay time - the time 
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FIG. 1: (Color online) Basic features of sacrificial bonds and 
hidden length, (a) Hypothesized structure of polymeric glue 
strands, with sacrificial bonds (blue circles) along the polymer 
chain backbone resisting chain rupture as they are stretched 
during microcracking. There are three types of sacrificial 
bonds: bonds within a polymer chain; bonds between dif- 
ferent chains; and bonds connecting the polymer chain to 
the substrate (collagen fibrils in the case of animal bone). 
Adapted from [4j. (b) Force change associated with sacrificial 
bond breakage and hidden length release, (i) Before a sacri- 
ficial bond is broken, only the black length of the molecule 
contributes to the entropic configurations and to the force 
with which the molecule resists stretching. The red length 
of the molecule is hidden from the force by the sacrificial 
bond, (ii) When the bond breakage threshold is reached, the 
bond breaks and the whole length of the polymer (black plus 
red) contributes to its configurational entropy. This sudden 
increase in entropy leads to an abrupt force drop, (iii) As 
the polymer molecule is further stretched, the force it sup- 
ports increases, until the entire molecule detaches from the 
substrate and ruptures. The grey area represents the extra 
work done in stretching a polymer with sacrificial bonds and 
hidden lengths, relative to a polymer of the same length but 
without such structural features. Reprinted with permission 
from 1141. 



that permits self-healing between two successive pulling 
experiments - on the macroscopic mechanical response 
of a pair of separating collagen fibrils. Based on these 
results, we propose in Sec. V a constitutive law that de- 
scribes the macroscopic response of separating collagen 
fibrils. 



II. KINETIC MODEL 

In this Section we introduce a kinetic model that de- 
scribes the dynamics of formation and breakage of sacrifi- 
cial bonds and release of hidden length in a polymer net- 
work, and relate this to mechanical forces on the polymer 
network. As a simplifying assumption, we propose that 
the network of polymers - in the specific example of an- 
imal bone, glue connection that hold the collagen fibrils 
together - can be described as a quasi-one-dimensional 
ensemble of polymer chains, regardless of whether sacri- 
ficial bonds are found within a single polymer molecule 
or between adjacent polymers. Each polymer chain may 
consist of one long polymer molecule or multiple poly- 
mers entangled together; see Fig. [TV a). The ensemble of 
these quasi-one-dimensional chains thus acquires a dis- 
tribution of total lengths. 

In addition, we assume that each polymer chain is 
semiflexible, so that the force-extension relation is given 
by the wormlike-chain model [T51 [15], in conformation 
with much of the literature on the mechanics of proteins: 



F = 



k B T 
b 



x 1 1 

L + 4(1 -x/L) 2 ~ 4 



(2.1) 



Here, the force F is entropic in nature, arising from the 
tendency of the polymer chain to recoil and return to a 



state of higher entropy as it is stretched. In Eq. (2.1), k B 
is the Boltzmann constant, T is the temperature, b is the 
persistence length, x is the end-to-end distance, and L is 
the contour length available for stretching - i.e. the to- 
tal contour length of each chain minus the hidden length 
shielded by sacrificial bonds. Breakage of each sacrificial 
bond unveils hidden length, resulting in a step jump in 
the available contour length L. This results in an in- 
crease in the chain entropy which causes abrupt force 
drops. Stretching the chain without breaking sacrificial 
bonds reduces the entropy, thereby dissipating a signifi- 
cant amount of energy. 

A sacrificial bond breaks when the force on the poly- 
mer chain exceeds the strength of that bond. As men- 
tioned in the Introduction, we assumed in our previ- 
ous work [14] that the bond strength is a uniform ran- 
dom variable, reflecting the randomness of bond break- 
age events. One approach to modeling the dependence 
of the mechanical behavior on pulling rate in the con- 
text of the previous model would be to represent the 
bond strength distribution itself directly as a function 
of pulling rate. However, this crude approach neglects 
the fundamental physics of rate dependence. Meanwhile, 
Bell's model [TS] expresses the transition probability of 
bond formation and breakage events as a Boltzmann fac- 
tor that involves the product of the force and a parameter 
with the dimensions of length, interpreted as the distance 
from the transition state of the conformational change. 
Along with other more sophisticated models (such as the 
Kramer theory [20]). it has been successful in account- 
ing for the rate dependence of forced protein unfolding, 
which in most cases involve breakage of weak internal 
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bonds. How, then, are we to apply such kinetic models 
to the breakage of sacrificial bonds and release of hidden 
length? 

We proceed by assuming that the breakage of sacrifi- 
cial bonds follows a two-state pathway, so that we can 
apply Bell's theory. At large forces and pulling rates, 
the formation of sacrificial bonds can be neglected. Mo- 
tivated by Su and Purohit [3T], we propose that the rate 
of change of the number of sacrificial bonds Nb is given 
by the first-order differential equation 



dt 



-k f N b 



(2.2) 



N£ is the continuous version of the integer Nb] it repre- 
sents the number of sacrificial bonds at a given instant of 
time, averaged over an ensemble of many polymer chains. 
It will be thresholded below (see Eq. (2.5 1 ) to isolate in- 



dividual bond breakage and formation events, and it co- 
incides with Nb whenever it is an integer. Nf = N — 2Nb 
is the number of free sites, with N — L/b being the num- 
ber of sites, or persistence lengths, in the polymer, kf 
and kb are the rates at which bond fragmentation and 
bond formation events occur; according to Bell's theory, 
they are given by 



k f 

h 



ao exp 



/3 exp 



f FAx 



V k B T 
FAx b \ 
k B T J 



(2.3) 
(2.4) 



F = F(x) is the force-extension relation given by 
Eq. ( |2.1| . Axj and Axb are the distances to the tran- 
sition state; ao and /3o are, respectively, inverse time 
scales which describe the rate at which bond breakage 
and formation events occur at zero pulling force. We 
have mentioned in the Introduction that the physics of 
sacrificial bond breakage and protein unfolding are sim- 
ilar. Based on parameter estimates for the unfolding of 
proteins in |21) . Axf, Axb and b are expected to be of 
the order of 0.1 nm. 

A bond breakage event occurs when Nb decreases by 
unity - that is, when N£ reaches an integer. Thus, the 
condition for a bond formation event to happen is 



dN* = / (-k f N b + k b N f )dt = 1 



(2.5) 



where the integral on the right hand side is over the time 
between successive bond formation events. Similarly, the 
condition for a bond breakage event to happen is 



J dNb' = J {-k f Nb + k b N f )dt 



= -1 



(2.6) 



where the integral on the right hand side is over the time 
between successive bond breakage events. In particular, 
for pulling experiments at constant velocity v, the pre- 
ceding equation gives 

/ 2 (k f (F(x)) N b - k b (F(x)) N f ) dx = v, (2.7) 



where x\ and x^ are the chain end-to-end distances at 
successive bond breakage events. 

In mechanical experiments on stretched glue connec- 
tion in animal bone [2r0] it has been found that sacri- 
ficial bonds mediated by ions such as calcium also form 
between the glue strand backbone and the collagen fib- 
rils. Breakage of these end bonds causes the detachment 
of the entire glue strand from one of the collagen fibrils, 
so that the stretching force on the glue strand imme- 
diately drops to zero. In addition, it has been found 
that broken links may self-heal, in that some broken end 
bonds would be restored if, after a particular pulling ex- 
periment, the entire sample is left untouched for times as 
short as a few seconds [4] . We propose that the breakage 
and restoration of end bonds can be described within the 
same theoretical framework. That is, the change in the 
number of end bonds N e is governed by the rate equation 



dt 



= -kf d N e 



+ kt nd (l-N e ). 



(2. 



N*, the continuous version of the integer N e , varies be- 
tween zero and unity; it can be interpreted as the fraction 
of end bonds that have yet to break. Notice that (1 — N e ) 
appears because each glue strand either attaches to the 



bone fibril, or does not anchor to it. 



K f 



and k? are 



the rates of end bond breakage and formation, necessarily 
different from kf and kb, given by 



uend 



a e exp 



P e exp 



' FAxf d " 
knT 



k R T 



(2.9) 
(2.10) 



As before, a e and /3 e are respectively the rates at which 
end bonds break and form when no external force is 



present, and Aa 



and Ax e b na 



are the distances from 



the transition state for end bond breakage and formation 
events. For a single polymer chain, the end bond breaks 
when 



kf d (F(x))dx = v 



(2.11) 



where x c is the chain end-to-end distance at which the 
end bond breaks and the chain detaches. 

We note on passing that for a collection of polymers 
stacked in parallel, Eq. (2.8 1 should more properly be 



interpreted as the governing equation for the fraction of 
polymer chains with end bonds restored as a function of 
time t: thus, 



dt 



kf Nl + kb (1 — N*) 



(2.12) 



Suppose t = marks the time at which all polymers 
detach from the surface, after the previous stretching ex- 
periment. Then the fraction A"* of polymers that adhere 
to the surface at time t is given by 



N* 



Pe 



1 



-(a e +/3 e )t 



(2.13) 



4 



Equations (2.7) and (2.11) can be used to predict the 



force-extension curve of a stretched polymer; in partic- 
ular, they can predict the chain end-to-end distance at 
which bond breakage events occur, as well as the corre- 
sponding bond strengths. Meanwhile, Eq. (2.13) is par- 



ticularly useful for analyzing the dependence of the me- 
chanical behavior of multiple polymers stacked in parallel 
on the delay time between pulls. 

Note that while our model is deterministic and does 
not capture the randomness of bond breakage events as 
in [14] . it represents an average over a large ensemble of 
experiments. 
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III. PULLING A SINGLE POLYMER CHAIN: 
THEORETICAL PREDICTIONS 

We begin by considering the force-extension behavior 
of a single polymer chain whose total length is L c 
1 00 nm. Let m denote the number of sacrificial bonds. 
For simplicity, we assume that shielded lengths do not 
contain sacrificial bonds, and that the length Lj of each 
hidden loop is a uniform random number less than L c /m. 
Then the initial available length is Li — L c — Ylj=i 
To locate bond breakage events, we integrate Eq. ( |2.7[ ) 
over the force-extension profile, and compute the end-to- 
end distance xi at which each individual bond breakage 
event occurs, assuming that bonds break in the order of 
increasing shielded length. We integrate Eq. (2.11) over 



the entire force-extension curve to locate the maximum 
pulling distance before the polymer chain detaches from 
the underlying material. 

Figure |2(a) show the force-extension curves of a poly- 
mer chain with m — 6 sacrificial bonds, stretched at four 
representative velocities v = 10 2 , 10 3 and 10 4 nm s _1 . In 
computing these theoretical curves we have used the pa- 
rameter estimates ao = 0.3 s _1 , 0q = 0.003 s™ 1 , a e = 0.1 
s _1 , b = 0.1 nm, Axf = 0.25 nm, Axb — 0.1 nm, and 
Ax e — 0.15 nm. The magnitudes of these parameters are 
roughly consistent with those in protein unfolding models 
(for example [21] )■ The average internal bond strength 
varies from roughly 80 pN at v = 10 2 nm s" 1 to 150 pN at 
v = 10 4 nm s~ 4 , while the end bond strength varies from 
roughly 150 pN to 300 pN over the same range of pulling 
velocities. For comparison, Fig. |2(b)] shows sample exper- 
imental pulling curves of a polymer chain with sacrificial 
bonds and hidden length. The force drops due to break- 
age of sacrificial bonds and release of hidden length are 
also seen in the grey curve here, which represents the be- 
havior of a typical polymer chain with sacrificial bonds, 
stretched by the tip of an atomic force microscope. The 
black curve shows the mechanical response of an oth- 
erwise identical polymer chain, i.e., one with the same 
total length but with no sacrificial bonds. In both fig- 
ures, the shaded area between the two curves represents 
the increase in toughness due to the presence of sacrificial 
bonds and hidden length. Due to the inherent variability 
of experimental parameters (ion concentration in buffer, 
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FIG. 2: (Color online) (a) Force-extension curves of a poly- 
mer with a total contour length of L c = 100 nm and with 
m — 6 sacrificial bonds, stretched at v = 10 2 nm s _1 (blue), 
10 3 nm s _1 (red), and 10 4 nm s _1 (green), increasing from 
bottom to top. (b) Comparison between sample experimen- 
tal pulling curves of a polymer chain with sacrificial bonds 
and hidden length (upper grey curve) and one without (lower 
black curve). Adapted from [2]. In both figures, the shaded 
area represents the extra energy dissipated through sacrificial 
bonds, a measure of the increase in toughness of the material. 



chain length, etc.) in different samples, the aim of this 
comparison with individual chain pulling experiments is 
to seek qualitative, rather than exact quantitative, agree- 
ment; qualitative agreement between the theoretical and 
experimental results is clear. 

Figure 3(a) shows 
stretch x 



the variation of the maximum 
as a function of the pulling velocity, for two 
polymer chains with 6 and 12 sacrificial bonds respec- 
tively but are otherwise identical; the overlap of the 
curves indicates that this is independent of the number of 
sacrificial bonds. Figure [3(b) | is a log-linear plot that dis- 
plays the variation of the end bond strength as a function 
of the pulling velocity v. The linearity of the plot shows 
that the end bond strength, which equals the maximum 
stretching force along the stretching profile, varies loga- 
rithmically with the stretching velocity v, as is predicted 
by Bell's theory [15] and seen in experiments with dentin 
matrix protein 1 [16] , shown here in Fig. |3(c) 
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FIG. 3: (Color online) (a) Maximum pulling distance, and 
(b) end bond strength, as functions of pulling velocity v, for 
m = 6 (blue circles) and 12 (red squares) sacrificial bonds. 
The overlap of the curves indicates that both quantities ap- 
pear to be independent of m. Figure (c) is the dynamic force 
spectrum of dentin matrix protein 1 with sacrificial bonds 
mediated by sodium, calcium and lanthanum buffers respec- 
tively, adapted with permission from [16j . Our theory pre- 
dicts a log-linear dependence, as observed in experiments, of 
the end bond strength on the pulling rate. 
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FIG. 4: (Color online) Total energy dissipation over the entire 
course of stretching, as a function of the number of sacrificial 
bonds m, for stretching velocities v — 10 2 nm s _1 (blue), 10 3 
nm s _1 (red), and 10 4 nm s" 1 (green), increasing from bottom 
to top. Results are averaged over 200 runs. The vertical bars 
indicate one standard deviation. 



Figure [4] shows the total energy dissipation, a mea- 
sure of the toughness of the polymer chain given by the 
area under the force-extension curve, as a function of 
the number of sacrificial bonds m, for three representa- 
tive stretching velocities v — 10 2 ,10 3 and 10 4 nm s _1 , 
and averaged over 200 runs. While sacrificial bonds and 
hidden length constitute a major toughening mechanism, 
increasing the number of sacrificial bonds beyond to 15 
fails to further stiffen the chain, as is found in [Hj. In 
addition, the importance of this toughening mechanism 
becomes more pronounced at high pulling velocities (of 
the order of 1000 nm s _1 ), providing increased resistance 
against impact loading. 



IV. STRETCHING MULTIPLE POLYMER 
CHAINS IN PARALLEL; EFFECT OF DELAY 
TIME BETWEEN PULLS 



In this section we consider the dynamical behavior 
of multiple polymers stretched in parallel. For simplic- 
ity, assume that each polymer chain is independent of 
the others, with no entanglement between the polymer 
strands, so that the total force equals the sum of forces 
in each polymer chain. 

Figure [5] shows the force-extension curves of N p — 200 
parallel polymer chains, at pulling velocity v — 1000 nm 
s , for delay times (the time between rupture of all end 
bonds and the start of the next pulling experiment) rang- 
ing from 1 to 20 seconds. In computing these theoretical 
curves we have used Eq. (2.13) to calculate the fraction 



of polymers with restored end bond connections to the 
underlying substrate (for example, mineralized collagen 
fibrils in the case of glue connection) as a function of the 
delay time t. We assume that the total contour length L c 
of each polymer is uniformly and randomly distributed 
between 50 and 150 nm. We use Eq. (2.5) to compute 
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FIG. 5: (Color online) Force-stretch curves for N p = 200 
polymer chains stacked in parallel. The delay times are t = 2 s 
(blue), 5 s (red), 10 s (green) and 20 s (orange) from bottom to 
top, corresponding N P N* = 8, 18, 28 and 36 polymer chains 
adhering to both pieces of substrate (collagen fibril) at the 
beginning (see Eq. ( |2.13[ ) for an expression for N*). 
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the number of internal sacrificial bonds that form as a 
function of the delay time, and assign random hidden 
lengths to each of these sacrificial bonds as before. As 
before, we have chosen a a = 0.3 s _1 , /3 = 0.003 s _1 , 
and a e = 0.1 s _1 . To account for the relatively slow re- 
covery of ruptured polymer chains as a function of the 
delay time, as seen in [4], we choose /3 e = 0.025 s _1 . 
These sample pulling curves indicate that the extension 
at maximum force and the maximum extension at a given 
pulling velocity is independent of the delay time. Also, 
the force p eaks level off for large delay times. This follows 
from Eq. (2.13), which shows that the fraction of poly- 
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mers with restored end bond connections approaches the 
asymptotic limit /3 e /(a e +/3 e ) as the delay time t becomes 
large; only those polymers with restored end bonds carry 
the pulling force and contribute to energy dissipation. In 
the limit a e -C j3 e , all end bonds are restored for large 
delay times between successive pulling experiments. Our 
present choice of parameters, however, stipulates that at 
most one-eighth of all glue strands are attached to bone 
fibrils at both ends. 



To verify these claims, Fig. 6(a) shows the peak force 
as a function of delay time t between pulls, for N p = 200 
parallel polymers pulled at a velocity v = 1000 nm s _1 , 
averaged over 100 runs. The peak force increases roughly 
in proportion to N* = (/3 e /(a e + (3 e ))(l - e ~(^+Mt) : 
the fraction of polymers that possess restored end con- 
nections and therefore transmit the force. Figure |6(b)| 
shows the displacement at maximum stretching force, as 
well as the maximum stretch, as a function of delay time 
between pulls. Both quantities are roughly independent 
of the delay time, except that the displacement at max- 
imum stretching force shows a slight decrease for small 
delay times. This can be traced to the fact that the small 
number of internal sacrificial bonds that restore for small 
delays times leads to fewer cusps in the force-extension 
curve of each polymer chain and reduces the extra energy 
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FIG. 6: (Color online) (a) Peak force, and (b) displacement at 
maximum force (blue, bottom curve) and maximum extension 
(red, top curve), as functions of delay time, for N p = 200 
polymers stacked in parallel, pulled at v = 1000 nm s _1 . 
Results are averaged over 100 runs, and the error bars indicate 
one standard deviation. 



dissipation brought by these sacrificial structures (see be- 
low), thereby delaying the occurrence of the maximum 
force. 

Figure [7(a) shows the total energy dissipation, a mea- 
sure of the toughness of the glue connection between 
the two pieces of underlying material, as a function of 
the delay time t between pulls. For small delay times, 
the marked growth in the number of restored end bond 
connections, the number of restored internal sacrificial 
bonds, and the increase in energy dissipation as a func- 
tion of the number of internal sacrificial bonds (cf. Fig. [4]) 
all contribute to the fast increase of total energy dissipa- 
tion as a function of the delay time. For large delay 
times, however, these growths slow down gradually, so 
that the increase in energy dissipation flattens out. A 
comparison to Fig. |7(b)| which shows the experimental 
measurements for the energy dissipation involved in the 
separation of two bone fibrils [3] , indicates that our theo- 
retical prediction qualitatively matches the experimental 
observations of bone fibrils in a buffer of calcium and 
sodium ions. The result for a sodium buffer corresponds 
to different choices for the rate parameters kf, kb, ctQ and 
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FIG. 7: (Color online) (a) Total energy dissipation over the 
entire course of stretching, as a function of delay time, for 
N p = 200 polymers stacked in parallel, pulled at v — 1000 nm 
s _1 . Results are averaged over 100 runs, and the error bars 
indicate one standard deviation, (b) Figure 2(c) from [4], re- 
produced here for comparison, shows the energy dissipation 
involved in the separation of bone fibrils, in a buffer where cal- 
cium and sodium ions are present (red) , and in a buffer where 
only sodium ions are present (blue). The authors there con- 
cluded that calcium ions lead to enhanced bond strength, and 
within our choice of parameters, our theoretical prediction for 
the energy dissipation qualitatively matches that of [4] in the 
presence of calcium ions. 
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V. TOWARDS A CONSTITUTIVE LAW FOR 
FORCE AS A FUNCTION OF DISTANCE 

In multiscale simulations of bone fracture it is neces- 
sary to incorporate a constitutive law for the separation 
of collagen fibrils under tensile stress. More specifically, 
in realistic situations where hundreds of glue strands are 
present between each pair of collagen fibrils, we need, in 
a mean-field sense (i.e. smoothing out all abrupt force 
drops due to bond breakage or detachment), a force law 



F(x,v,t) for the force on the polymeric system under 
stretch. 
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FIG. 8: (Color online) Force per polymer chain that is at- 
tached at both ends at the beginning of the experiments, as a 
function of pulling distance, for N p = 200 polymers. Results 
are averaged over 100 runs, thereby smoothing out abrupt 
force drops. The delay times are t — 5 s (blue), 10 s (red) and 
20 s (green). The pulling rates are v = 100, 1000 and 10000 
nm s _1 for the three sets of curves, from bottom to top. 

To this end, Fig.[8]shows sample force-extension curves 
of N p = 200 polymer chains stacked in parallel, normal- 
ized by the number of chains that are attached at both 
ends at the beginning, ignoring all interactions between 
them, for delay times ranging between 1 and 10 seconds. 
In computing these theoretical curves we have averaged 
over 100 pulling experiments and divided the total force 
F(x,v,t) by the total number of polymer chains N p N* 
that are attached to the collagen fibrils at both ends at 
the beginning, where N* = {j3 e /{a e +l3 e )){l-e~ { - a ' +l3 ^ t ) 



as given by Eq. (2.13|. As in Fig. [5] the total contour 
length L c of each polymer is uniformly distributed be- 
tween 50 and 150 nm. Importantly, for each pulling ve- 
locity v, the curves for different delay times t collapse to- 
gether. This implies that in the limit of long delay times 
t, the total force in separating two pieces of collagen fib- 
ril is given as a function of distance x and separation 
velocity v by 

F(x,v,t) = N p N* e (t)f(x,v). (5.1) 

The delay time dependence comes in only through the 
fraction N*(t) of polymer chains attached to the bone 
fibrils at the beginning. The average force on each of 
these polymer chains, f(x,v), is independent of the de- 
lay time t and can be approximated by separate power 
law fits to the increasing (strengthening) and decreasing 
(weakening) portions of the curves: 



f(x,v) 



f P (v) 

f P (v) 
10, 



x p (v) 
x c (v) 



x c (v) - x p (v) 



if x < x p (v); 

if x p (v) < x < x c (v); 
if x > x c (v). 

(5.2) 



Here, f P (v), x p (v) and x c (v) are the velocity-dependent 
peak force, end-to-end distance at peak force, and max- 
imum pulling distance, respectively. We find that s x » 
1.35 and S2 ~ 0.65 and, in the case v — 1000 nm s _1 
shown in Fig. [Sj that f p 0.094 nN, x p 43.5 nm and 
x c 115 nm. Figure [9] shows how this function fits sam- 
ple force-displacement profiles. 




FIG. 
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9: (Color online) Force-stretch curves as in Fig. [H] for 
N p = 200 polymers stacked in parallel, pulled at v = 1000 
nm s , normalized by the number of polymer chains that are 
attached at the beginning of the experiment. Fhe delay times 
are t — 2 s (blue), 5 s (red), 10 s (green) and 20 s (orange). 
The thick black curve represents the fitting function Eq. ( |5.2[ ) 
with parameter values given in the text. Note that the blue 
curve for the short delay time t — 2 s is noticeably jagged, 
owing to the fact that only 8 polymer chains (see Eq. (2.13 1 ) 
adhere to the bone fibrils at the beginning. The quantities 
f p , x p and x c appearing in Eq. (5.2| are labeled in the figure. 



The quantity x„(v) marks the transition from a 
strengthening to a weakening behavior, associated with 
the gradual detachment of polymer chains. To check this 
assertion, the solid curves in Fig. [lO] show the evolution 
of the fraction of initially intact glue strands that remain 
attached to the collagen fibrils at both ends (i.e. with end 
bonds being intact), as a function of the displacement x, 
at the same pulling velocity v — 1000 nm s _1 . One sees 
immediately that the onset of polymer chain detachment 
coincides with the transition to weakening behavior at 
x p . In addition, under our assumption of a uniform dis- 
tribution of polymer chain lengths, the number of poly- 
mer chains that have yet to rupture decreases linearly 
with separation distance in the weakening regime. On 
the other hand, the dashed curves show that breakage 
of sacrificial bonds within individual polymer chains oc- 
curs continually in both the strengthening and weakening 
regimes. Figure 10 thus demonstrates how microscopic 



physics (bond breakage) , characterized by the state vari- 
ables Nf, and N e , accounts for macroscopic behavior (dis- 
placement strengthening and subsequent weakening) in 
biological structures that contain sacrificial bonds and 
hidden length. In the context of dynamic fracture, the 
strengthening regime corresponds to crack arrest and the 
weakening regime corresponds to crack propagation and 



catastrophic failure. 
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FIG. 10: (Color online) Number of polymer chains that have 
yet to rupture, for N p = 200 polymers stacked in parallel, 
pulled at v = 1000 nm s" 1 , normalized by the number of 
polymer chains attached to bone fibrils at the beginning (solid 
curves); and fraction of internal sacrificial bonds that have 
yet to rupture (dashed curves). The delay times are t = 2 s 
(blue), 5 s (red), 10 s (green) and 20 s (orange). The black 
dot-dashed curve is a vertically rescaled f(x, v) and shows 
that the onset of polymer chain detachment coincides with 
the transition to weakening behavior after reaching the max- 
imum force fp. This illustrates the interplay between mi- 
croscopic physics (bond breakage) and macroscopic behavior 
(displacement weakening). 

To extract the velocity dependence of the quantities 
f P (v), x p (v) and x c (v) explic itly, we plot these quanti- 
ties in Figs. |ll(a) and 11(b) They can be fit with the 
functional forms 



f P (v) 
x p (v) 

x c (v) 



/ilog( - ] +/o; 



Pi log 



Po; 



(■■2 



log ( — 



Ci log — 
>0 



(5.3) 
(5.4) 

c . (5.5) 



100 

41 nm, 
= -0.57 



Within our choice of system parameters, we find Do 
nm s~\ f = 0.064 nN, f x = 0.013 nN, p = 
Pi = 1.3 nm, Co = 104 nm, C\ — 6.3 nm and C2 
nm. 

We have thus shown that in the limit of long delay 
times t, the total force F on the ensemble of polymers 
can be factored into the product of the number N p N* (t) 
of polymer chains that are intact at the beginning of the 
pulling experiment, times the average force f{x,v) per 
polymer chain. Both the strengthening and weakening 
regimes - the latter being associated with the rupture of 
polymer chains and their detachment from the substrate 
- can be described by power laws characterized by the 
velocity-dependent peak force f p {v), the displacement 
at peak force x p (v) and the maximum pulling distance 
x c (v). Among these, the peak force f p varies linearly 
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FIG. 11: (Color online) (a) Plot of the average peak force 
f p (v) per polymer chain as defined in Eq. ( |5.2[ ) versus the 
pulling speed v. The blue data points are average values 
from 10 ru ns and the red curve is the log-linear fit given 
by Eq. (5.31. (b) Plots of the displacement at peak force 



x p (v) (blue circles) and maximum pulling distance x c (v) (red 
squares) as defined in Eq. ( |5.2[ ) versus the pulling speed v. 
The data points are average values from 100 runs. The black 
curves are the fits given by Eqs. (5.4 1 and (5.51. 



with the logarithm of the pulling velocity v, in confor- 
mity with experiments. This constitutive approach en- 
ables us to describe the mean-field dynamics of sacrificial 
bond breakage and hidden length release using several 
adjustable parameters, without having to account for the 
random breakage of individual bonds in detail. 



Nb and the number of chains N e that adhere to the 
substrate - the only two molecular-state-dependent vari- 
ables in the theory - with the pulling distance x and 
the total force F on the polymer network, according to 
a velocity-dependent criterion that determines the times 
or displacements at which bond breakage occurs. The 
force is entropic in nature, given by the wormlike-chain 
model as a function of the pulling distance x and the 
amount of available contour length L, the latter of which 
is computed in terms of the number of remaining sacrifi- 
cial bonds Nf,. 

We have shown that sacrificial bonds and hidden 
length lead to a marked increase in fracture toughness 
in materials where they are present, and that both the 
fracture toughness and maximum displacement before 
complete rupture in a pulling experiment increases with 
the pulling velocity v which drives the system away from 
equilibrium. In particular, the peak force f p in the force- 
displacement profile varies linearly with the logarithm of 
the pulling velocity v, in conformity with various mechan- 
ical experiments on biological molecules such as [16] . In 
addition, our kinetic model naturally incorporates self- 
healing, evidenced by the increase in the number of at- 
tached polymer chains, rupture peak height and total 
fracture toughness with recovery time. Our simple quasi- 
one-dimensional model, however, does not explicitly ac- 
count for the effect of crosslinks and entanglements in a 
network of glue strands. The extent to which these ad- 
ditional microscopic mechanisms will impact the macro- 
scopic behavior will be investigated in future work. 

Based on our theoretical calculations we have pro- 
posed a phenomenological description for the force- 
displacement profile of a collection of polymer chains 
with a distribution of lengths. The force-displacement 
profile consists of a strengthening regime for small dis- 
placements, where the force increases with the displace- 
ment according to a power law. This is followed by a 
weakening regime associated with the gradual detach- 
ment of polymer chains which no longer contribute to 
force transmission. Such a constitutive description will 
be of utmost utility in future multiscale simulations of 
bone fracture. The dynamical behavior of glue connec- 
tion between collagen fibrils has important implications 
on crack propagation, crack arrest, strength recuperation 
and collagenous diseases in bone. 



VI. CONCLUDING REMARKS 
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